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Abstract — We show how good quantum error-correcting codes 
can be constructed using generalized concatenation. The inner 
codes are quantum codes, the outer codes can be linear or 
nonlinear classical codes. Many new good codes are found, 
including both stabilizer codes as well as so-called nonadditive 
codes. 

Index Terms — Generalized concatenated codes, quantum error 
correction, stabilizer codes, nonadditive codes 

I. Introduction 

The idea of concatenated codes, originally described by 
Forney in a seminal book in 1966 [10], was introduced to 
quantum computation three decades later [1], [12], [18], [19]. 
These concatenated quantum codes play a central role in fault 
tolerant quantum computation (FTQC) as well as in the study 
of constructing good degenerate quantum codes. 

Blokh and Zyablov [3], followed by Zinoviev [25] intro- 
duced the concept of generalized concatenated codes. These 
codes improve the parameters of conventional concatenated 
codes for short block lengths [25] as well as their asymptotic 
performance [4]. Many good classical codes, linear and non- 
linear, can be constructed using this method. 

In [17] we, together with Smith and Smolin, have introduced 
generalized concatenated quantum codes (GCQC). It is shown 
that GCQC in its simplest form, i. e., two level concatenation, 
is already a powerful tool to produce good nonadditive quan- 
tum codes which outperform any stabilizer codes. 

This paper focuses on the multilevel concatenation for 
quantum codes. We use the framework of stabilizer codes and 
the generalization to codeword stabilized (CWS) codes [6], 
[7] and union stabilizer codes [15], [16]. This allows to use 
classical codes as outer codes. We further extend our multilevel 
concatenation technique to the case of different inner codes, 
which allows us to construct codes of various lengths. 

II. Background and Notations 

A general quantum error-correcting code (QECC), denoted 
by C = ((n, K, d)) q , is a K -dimension subspace of the Hilbert 
space Hf n of dimension q n that is the tensor product of n 



complex Hilbert spaces H q = C q of dimension q. Here we 
restrict q = p m to be a prime power. A QECC with minimum 
distance d allows to correct arbitrary errors that affect at most 
(d — l)/2 of the n subsystems. 

Most of the known QECCs are so-called stabilizer codes 
introduced independently by Gottesman [11] and Calderbank 
et al. [5]. The code is defined as the joint eigenspace of a set 
of commuting operators [11]. Equivalently, the code can be 
described by a classical additive code C over GF(q 2 ) that is 
self-orthogonal with respect to a symplectic inner product [2], 
[5]. Denoting the symplectic dual code by C*, the minimum 
distance of the quantum code is given by 

d = min{wgt(c) :ceC*\C}> d min (C*). 

If d = G? m ; n (C*), the quantum code is called pure or non- 
degenerate. The corresponding stabilizer (or additive) code is 
denoted by C = [[n, k, d]] q and has dimension K = q k . 

The first nonadditive code ((5, 6, 2))2 which has a higher 
dimension than any stabilizer code of the same length correct- 
ing one erasure can be explained as the union of six locally 
transformed copies of the stabilizer code [[5,0, 3]]2 (see [14], 
[22]). A one-dimensional stabilizer code [[n, 0, d]] can also be 
described by a graph with n vertices [23]. The corresponding 
quantum states are referred to as graph states. Combining 
locally equivalent graph states, the first one-error-correcting 
nonadditive quantum code ((9, 12, 3))2 with higher dimension 
than any stabilizer code has been found [24]. The theoretical 
ground for these codeword stabilized (CWS) quantum codes 
has been laid in [6], [7]. 

In [15], [16], the framework of union stabilizer codes 
has been introduced. Starting with a stabilizer code Co = 
[[n, k, do]] q , a union stabilizer code is given by 

C = tC , 

t£T 

where To = {ti, . . . , tx} is a set of tensor products of (gener- 
alized) Pauli matrices such that the spaces tiCo are mutually 
orthogonal. Then the dimension of the union stabilizer code 
C is Kq k , and we will use the notation C = ((n, Kq k ,d)) q . 
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Similar to stabilizer codes, a union stabilizer code can be 
described in terms of classical codes. Given the symplectic 
dual Cq of the additive code Co associated to the stabilizer 
code Co, the union normalizer code is the union of cosets of 
Cq given by 

C*= IJ (%+t = {c + t j :ce(%,j = l,...,K}. (1) 

teT 

Here T is the set of vectors tj G F^ 2 corresponding to the 
generalized Pauli matrices f$ G Tq. 

Proposition 1 (cf. [16]): The minimum distance of a union 
stabilizer code with union normalizer code C* is given by 

d = min{wgt(w) : v e (C* - C*) \ C } 
> d min (C*) 

= min{dist(c + U, c + U,ti> G Tq, c, c G Cq 

c + U^ c' + ti,}, 

where C* - C* := {a - b: a,b G C*} denotes the set of all 
differences of vectors in C*, and Co < Cq is the symplectic 
dual of the additive closure of the (in general nonadditive) 
union normalizer code C*. 

Hence in order to construct a union stabilizer code with 
distance d, it suffices to find a large classical code C* with 
minimum distance d that can be decomposed into cosets of 
an additive code Cq that contains its symplectic dual. Two 
extremal cases are stabilizer codes where only one coset is 
used, and CWS codes for which Cq = Co is a symplectic 
self-dual code. 



the Hilbert space TCf N . Denote the standard basis of each 



H%» by 



{{if 1 ) eg) ... ® < »W < qj - 1, 1 < v < N} 

(where j runs from 1 to r), and the bases of the codes Aj are 
denoted by {\(f)f. ) : < lj < Mj — 1}. Expanding the basis 
vectors of Aj with respect to the standard basis of Ti.f N we 
obtain 
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The basis vectors of the tensor product of all outer codes are 
given by 



where lj runs from to Mj — 1. Expanding these basis vectors 
with respect to the standard bases we obtain 
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III. Generalized Concatenation 

The basic idea of generalized concatenated quantum codes 
[17] uses just two levels of concatenation. Here we first present 
multilevel concatenation for quantum codes. Then we discuss 
a special case that can be described by classical codes only. 

A. Multilevel Concatenation for Quantum Codes 

The inner quantum code B^ = ((n, q^q 2 ■ ■ ■ q r , d\)) q is 
first partitioned into qi mutually orthogonal subcodes 
(0 < ii < q% — 1), where each B^' is an ((n, q% ■ ■ ■ q r , d 2 )) q 
code. Then each B^ is partitioned into q 2 mutually or- 

(2) (2) 

thogonal subcodes B- ' (0 < i 2 < 92 — 1), where B- ' 
has parameters {{n,qs ■ • ■ q r ,ds))q, and so on. Finally, each 
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(r-2) 
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is partitioned into q T -\ mutually orthogonal sub- 



codes B 
Thus 
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((n,q r ,d r )) q for < i r _i < q r -i — 1. 
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fl(0) = ® fl « *£> = ©*&, 



(2) 



ii=0 



i 2 =0 



and d\ < d 2 < . . . < d r . A typical basis vector of B^ will be 
denoted by \ip ili2 ... ir ) (0 < ii < q\ - 1, . . . , < i r < q r - 1), 
with subscripts chosen such that \<Pi 1 i 3 ,,,i r ) is a basis vector 



of all BP , Bf) 



(r-l) 
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In addition, we take as outer codes a collection of r quantum 
codes A\, . . . , A r , where Aj is an ((N, Mj, Sj)) qj code over 
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The basis of the resulting generalized concatenated quantum 
code Q is given by replacing the basis vectors in Eq. © using 
the mapping 

for 1 < < AT. Hence the basis of Q is given by 

IVW.Jr) = E a l!W) l < 1 ) 



„(r) 



® 1^,(1), (2) f w>. 



So Q is a quantum code in the Hilbert space Hf Nn of 
dimension M = M\M 2 ■ ■ ■ M r . As already mentioned, the 
construction given in [17] is a two-level construction with 
r = 2, while the concatenation of quantum codes used in the 
context of fault tolerant quantum computation (cf. [1], [12], 
[18], [19]) is a one-level construction, i.e. r = 1. 

B. Classical Outer Codes 

From now on we restrict ourselves in constructing union sta- 
bilizer codes. For simplicity we consider only nondegenerate 
codes here. 
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We take the inner code B^°' to be an ((n, K q k , d\)) q nonde- 
generate union stabilizer code, given by a classical symplectic 
self-orthogonal additive code Co C Cq — (n,q n+k ,d r ) q 2 and 
a set 7"(°) of K = qxq 2 ■ ■ ■ q r -i coset representatives. The 
corresponding classical union normalizer code is 



C* = 23* (0) = [j C* + t. 



The decomposition d2J of the inner quantum code into 
mutually orthogonal union stabilizer codes is based on the 
decomposition of the union normalizer code B*' ' that is 
obtained by partitioning the coset representatives 
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This defines union normalizer codes <B*W given by 



B 



ZH 2 ■ ■ -lj -1 
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C n * + t. 



The coset representatives in T^°> will be denoted by ij 1 i 2 ...i r _ 1 
with < ii < qi — 1, < i r _j < g r _i — 1. 
The indices are chosen such that i^, j r _ 1 belongs to all 

-r(i) -r(2) -r^- 2 ) 

il ' il»2 ' ' ' ' ' *l«2...«i — 2 ' 

Here £>*(°' is a classical code over GF(q 2 ) with parameters 
(n,qiq2 ■ ■ ■ q r -iq n+k ,di) q 2 that is the union of qi disjoint 



codes B, 



*(i) 



k(r-2) 



nally, each S^..., 



(n, q 2 

2) 



B 



ir-l) 



Qr-iq n+k , cfe) g 2 , and so on. Fi- 
is the union of q r -i disjoint codes 
(n, q n+k , d r ) q 2, each of which is a single coset 



of the additive code Cq. 

In total we use r classical outer codes. For the first r — 1 
outer codes we take At — (N, Mi, Si) qi , a classical code over 
an alphabet of size q^ with length N, size Mi, and distance 
Si. The code A r is a trivial code A r — [N, N, l] qr where 
q r = |C *| = q n + k . 

Next we show how to construct the classical generalized 
concatenated code using the inner code and the outer 

codes A\, . . . ,A r - What follows is an adaption of [20, Ch. 18, 
§8.2]. The trivial classical code A r = [N, N, l] Qr on level r is 
concatenated with the additive normalizer code Cq, resulting 
in the additive code (Cq) n which contains its symplectic 
dual Cq^. Note that this corresponds to concatenating a trivial 
quantum code A r = [[N, N, l]]„k with the stabilizer code 
Co- As a technicality we note that the alphabet size of the 
trivial classical outer code A r is q n+k , while the trivial outer 
quantum code A r is over quantum systems of dimension q k . 

The first r — 1 outer codes are used to define a set of coset 
representatives. For this, form an N x (r — 1) array 



(1) (2) 
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N 
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where the first column is a codeword of Ai, the second is 
in Ai, etc. Then replace each row 



,(D (2) 

,u 3 



a) by 



the coset representative tji) n (2) (r- ij 



= Tj. (For this, 



J 7 J 

label the elements of the alphabet of size q,; by the numbers 
0, 1, . . . , qi — 1 in some arbitrary, but fixed way.) The resulting 
N x n arrays T = (Ti,...,Tjy) (considered as vectors of 
length Nn) form the new set of coset representatives of the 
generalized concatenated code 

C c = |J (C *x...xC *) + (Ti,...,T iV ). (5) 

(T U ...,T N ) 

Clearly, this code C* c has the form of a union normalizer code 
as specified in ([TJ. Hence C* c defines a QECC. The properties 
of this code are given by the following theorem. 

Theorem 2: The minimum distance of the union normalizer 
code 

C gc = {{nN, MtM 2 ■ ■ ■ M r ^ iq kN , d)) q , 

corresponding to C* c given in <(5j is 

d > min{(5i<ii, . . . , 8 r -id r -\, d r }. 

Proof: Let c and c be two distinct codewords of C* c . If 
they belong to the same coset, then c — c S (Cq) . Hence 
their distance is at least d r . Now assume that c and c lie 
in different cosets given by the arrays (a*- ) and (a,j ). If 
the arrays differ in the V th column then they differ in at 
least 8 V places in the v th column. By definition t^i 2 ,.,i 
and ii 1 < 2 ...i 1/ _ 1( 8... (with a =/= (3) both belong to 

Therefore the corresponding codewords of ( differ 

in at least d v places. Hence c and c differ in at least 5 v d v 
places. ■ 



(!/-!) 



C. Additivity Properties 

We know that if C* c is an additive code, then the correspond- 
ing quantum code C gc is a stabilizer code. So the question is 
when does generalized concatenation yield an additive code. 
The following is an adaption of a result from [9]. 

Proposition 3: Given additive, i. e., F p -linear, outer codes 
A\, . . . , A r —x and an additive inner code B, the resulting 
generalized concatenated code is additive if the mapping 



I (1) (2) 



Ar-1) 



) i ► t 



(6) 



is Fp-linear. 

Hence we can construct stabilizer codes from a sequence of 
nested stabilizer codes yielding a decomposition of the inner 
code and classical linear outer codes. 
Theorem 4: Let 

= [[n, ko, di)] q D S« = [[n, k u d 2 ]] q D . . . 

..oB'^ = [[n,kr-i,dr]] q 

be a sequence of nested nondegenerate stabilizer codes. This 
defines a decomposition of the inner code B^°\ Using r — 1 
additive outer codes Ai — (N,Mi,6i) qi where qi — q k i-~>-~ k i 
together with the trivial code A r = [N,N, l] qr where q r = 



n-\-k r —i 



, by generalized concatenation we obtain a stabilizer 



code with parameters [[nN, K, d]] q where 

d > min{(5iG?i, 5 2 d2, . . . , <5 r _id,-_i, d r } 
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and 

K = k?\og q (M 1 M 2 ---M r _ 1 ). 
Examples for this theorem are given in the next section. 

IV. Examples 

A. Stabilizer Codes 

Example 5: Consider the following sequence of nested sta- 
bilizer codes: 

B (0) = [[6, 6, l]] a D B« = [[6, 4, 2]] 2 D B^ = [[6, 0, 4]] 2 . 

The largest code B^ can be decomposed into four mutually 
orthogonal subspaces, each of which is a code [[6,4, 2]] 2 . 
Then each of these codes B^ is decomposed into 16 one- 
dimensional spaces [[6, 0, 4]] 2 . Hence we need nontrivial outer 
codes with alphabet sizes 4 and 16, which we chose to be 

Ai = [6,3,4] 4 and A 2 = [6,5,2]i 8 , 

together with A3 = [6, 6, l] 2 e. The dimension of the resulting 
code is I .Ai| x |A 2 | = 4 3 16 5 = 2 6 2 20 = 2 26 , and the minimum 
distance is at least min{4 x 1, 2 x 2, 4} = 4. Taking an additive 
map ©, we obtain a stabilizer code. As all inner codes 
are GF(4)-linear, we can even chose the mapping © to be 
GF(4)-linear, resulting in a GF(4)-linear code [[36, 26, 4]] 2 . 
This code improves the lower bound on the minimum distance 
of a stabilizer code [[36, 26, d]] 2 given in [13]. 

Our construction allows to adopt most of the known varia- 
tions of generalized concatenation for classical codes. In [8] a 
modified generalized concatenation has been introduced which 
uses outer code A, of different lengths m as well as different 
inner codes Bj°\ 

Example 6: Using the stabilizer code B^ = [[21, 15, 3]] 2 , 
we can decompose the full space B^ = [[21, 21, 1]] 2 into 64 
mutually orthogonal codes [[21, 15, 3]] 2 . In order to construct 
a generalized concatenated quantum code of distance three, 
we need a classical distance-three code over an alphabet of 
size 64, e.g., the classical MDS code Ai = [65,63, 3] 2 e, 
as well as the trivial code A 2 — [65, 65, l] 2 2i+is. Then by 
generalized concatenation one obtains a perfect quantum code 
[[1365, 1353, 3]] 2 . Instead of taking 65 copies of the inner 
code of length 21, we can use any combination of inner codes 
Bf ] = [[nj,nj - 6, 3]] 2 with nj £ {7, ... , 17, 21}. Note that 
now the trivial outer code A 2 has to be modified in such a 
way that by concatenation we get the normalizer code of the 
direct product of the various inner codes Bj \ Overall we 
obtain quantum codes with parameters [[n,n — 12,3]] 2 for 
n = 455, . . . , 1361 and n = 1365. 

Note that for quantum codes, the existence of a code 
[[n, k, d]] q does not necessarily imply the existence of a 
shortened code [[n — s,k — s,d\] q . In general, one would have 
to analyze the weight structure of an auxiliary code, the so- 
called puncture code, introduced in [21]. Varying the length 
of the inner quantum codes, we can directly construct shorter 
codes. 



B. Nonadditive Codes 

In our construction, we can also use classical nonlinear 
codes as outer codes. Good nonlinear codes can be obtained 
as subcodes of a linear code over a larger alphabet (or one of 
its cosets) by taking only those codewords whose symbols are 
taken from a subset of the alphabet. The following result can 
be found in [9, Lemma 3.1]): 

Proposition 7: If there exists an (n, K, d) q code, then for 
any s < q, there exists an (n',K',d) s code with size at least 
K(s/q) n . 

Example 8 (cf. [17]): We start with the sequence of inner 
codes 

B(°) = [[5,5,1]] 2 D5( 1 ) = [[5,1,3]] 2 . 

For the nontrivial outer code we take a code over an alphabet 
of size 2 5_1 = 16 and distance three. From the linear MDS 
code [18,16, 3] 17 over GF{11) we can derive a nonlinear 
code Ai = (18, [-^-],3)i 6 over GF(16) using Proposition 
Q The resulting generalized concatenated quantum code has 
parameters ((90, 2 81 825 , 3)) 2 , while the best stabilizer code 
has parameters [[90, 81, 3]] 2 . 

In the final example, we use three levels of concatenation and 
a nonlinear classical outer code. 

Example 9: Decompose the code = [[8, 8,1]] 2 using 
the sequence of nested stabilizer codes 

B (0) = [[8,8,1]] 2 D B« = [[8,6,2]] 2 D B^ = [[8,3,3]] 2 . 

As outer codes we need a code with alphabet size 2 8 ~ 6 = 4 
and distance three, a code with alphabet size 2 6-3 = 8 and 
distance two, as well as a trivial code. We take the nonlinear 
code Ai = (6, |~4 6 /5 2 ] , 3) 4 derived from the linear MDS code 
[6,4,3]5 over GF(5), the linear code A 2 = [6, 5,2] 8 over 
GF(8), and the linear code A3 = [6, 6, l] 2 s+3. The dimension 
of the generalized concatenated quantum code is |Ai | x |A 2 | x 
dim(_B^ 2 ^) 6 = 164 x 8 5 x 2 3x6 . Hence we get a nonadditive 
code ((48, 2 40 356 , 3)) 2 , which has a higher dimension than 
the best possible additive code [[48, 40, 3]] 2 . 

V. Decoding 

One of the advantages of concatenated codes as well as 
generalized concatenated codes is that decoding can be based 
on decoding algorithms for the constituent codes [9], [10]. For 
quantum codes, however, it is not possible to directly measure 
the "code symbols". Instead, decoding is based on measuring 
an error syndrome. 

For stabilizer codes, the error syndrome is obtained by 
measuring the eigenvalues of generators of the stabilizer group. 
The error syndrome can be defined in such a way that it 
corresponds to the error syndrome of the underlying classical 
code, and hence a classical decoding algorithm can be used. 

For generalized concatenated quantum codes derived from 
a sequence of nested stabilizer codes as in Theorem |4] the 
corresponding stabilizer groups are nested as well, with the 
stabilizer of the smallest code B^ being the largest. It is 
possible to choose its generators in such a way that stabilizers 
of the larger codes are generated by appropriate subsets. Hence 
the components of the syndrome vector reflect the nested 
structure of the inner code. 
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Again, we may not directly measure the syndromes of 
the N copies of the inner code. Instead, we compute the 
eigenvalues using some auxiliary quantum systems. Then we 
derive syndromes for the outer codes which will be measured. 

Details of the quantum circuits for syndrome measurement 
and iterative decoding algorithms are left to further work. 
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